Two-band superconductivity in doped SrTiC>3 films and interfaces 
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We investigate the possibility of multi-band superconductivity in SrTiOs films and interfaces us- 
ing a two-dimensional two-band model. In the undoped compound, one of the bands is occupied 
whereas the other is empty. As the chemical potential shifts due to doping by negative charge 
carriers or application of an electric field, the second band becomes occupied, giving rise to a strong 
enhancement of the transition temperature and a sharp feature in the gap functions, which is mani- 
fested in the local density of states spectrum. By comparing our results with tunneling experiments 
in Nb-doped SrTiOs, we find that intra-band pairing dominates over inter-band pairing, unlike 
other known multi-band superconductors. Given the similarities with the value of the transition 
temperature and with the band structure of LaAlOs /SrTiOs heterostructures, we speculate that the 
superconductivity observed in SrTiOs interfaces may be similar in nature to that of bulk SrTiOs, 
involving multiple bands with distinct electronic occupations. 
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The nature of multi-band superconductivity 
been the subject of intense experime 
theoretical (343 debate for decades. In the mid-1960s, 
specific heat measurements found evidence of multiple su- 
perconducting gaps in some elemental superconductors, 
such as Nb and V. More recently, multi-band supercon- 
ductivity has been found in compounds displaying rela- 
tively high transition temperatures T c : the magnesium 
diborides, with T c « 39K[Ioj], and the iron pnictides, 
with T c up to 56K pj] - [l6| . While the superconducting 
states in these two classes of materials have significantly 
different microscopic structures, many of their thermody- 
namic properties display common signatures associated 
with the existence of multiple gaps. Such features en- 
compass the temperature dependence of the specific heat 
and of the penetration depth, as well as the multiple-peak 
structure in tunneling spectroscopy. 
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Figure 1: (Color Online) A general illustration of the 
two band model as applied to bulk STO (left panel) and 
LAO/STO (right panel). For STO, the d xz / yz (band 1) crosses 
the Fermi level and the bottom of the d xy (band 2) is located 
above the Fermi level by \i c — /jl. In the case of LAO/STO 
interface, orbital reconstruction switches the position of the 
bands. 



Some of these characteristic multi-gap features have 
been observed in a lesser investigated material, electron- 
doped SrTiOs (STO). Indeed, oxygen- deficient STO pos- 
sesses a superconducting state below T c « 0.3 — 0.4 
K[l|, ISjl. Remarkably, tunneling measurements by Binnig 
et al. in Nb-doped STO in the early 1980s found two 
peaks in the local density of states beyond a certain elec- 
tronic density, providing strong evidence for multi-band 
superconductivity. More interestingly, the appearance of 
a second superconducting gap is correlated with a sudden 
increase in T c , which can reach values up to 0.7K. Band 
structure calculations on Nb-doped STO also finds that 
near the electronic concentration where T c is maximum, 
an additional electron-like band crosses the Fermi level 

Besides bulk STO, the recently discovered heterostruc- 
tures of LaA10 3 /SrTi0 3 (LAO /STO) also display super- 
conductivity up to T c ~ 0.2 K|2lL [22l|, which can be en- 
hanced to 0.3 — 0.4 K by the application of an electric 
field [23|, IU • This is surprisingly close to the aforemen- 
tioned T c value in bulk STO. Recent experiments have 
also found strong evidence that the maximum T c in the 
heterostructures is achieved once an extra electron band 
becomes occupied [17j|. In view of the early experiments 
by Binnig et al. [1], it is tempting to trace an anal- 
ogy between superconductivity in bulk STO and in the 
LAO/STO heterostructures. Thus, the understanding of 
the multi-band superconducting state in electron-doped 
STO may be relevant for the ongoing debate on the na- 
ture of the superconducting LAO/STO heterostructures. 

In both systems - bulk STO and LAO/STO interfaces 
- the origin and mechanism of the SC state is still an 
open issue. Most of the studies, however, have employed 
a single-band approach fid I25l ~ l28l |. In this paper, we do 
not discuss the origin of the SC state. We instead use a 
mean-field phenomenological two-dimensional model to 
focus on the consequences and signatures of multi-band 
superconductivity in these systems for varying electronic 
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concentration. Therefore, our approach is most relevant 
for thin films of doped STO and interfaces of STO with 
other systems that may provide negative charge carriers, 
such as LaA103 and interfacial oxide gel(2Q|. 

In our model, motivated by first-principle calculations 
d, the d xz ^ yz and d xy orbitals form two electron 
bands such that, in the weakly doped compound, only 
one of them crosses the Fermi level (see Fig. [1]). As 
the system is doped with negative charge carriers and 
the chemical potential \i shifts upwards, the second band 
crosses the Fermi level, causing important changes in the 
SC properties of the material. We calculate T c , the gap 
functions, and the local density of states (LDOS) for dif- 
ferent electronic concentrations, and compare our results 
with the tunneling experiments by Binnig et al. in Nb- 
doped STO Such a comparison demonstrates that an 
attractive intra-band (intra- orbital) pairing interaction 
dominates over the inter-band (inter- orbital) pairing, in 
contrast to other recently found multi-band superconduc- 
tors, such as the iron pnictides. 

Although small, the inter-band coupling causes a sharp 
increase of T c when the second band is occupied at 
fi = fjb c , and also enhances both gap functions, with 
the gap in the initially unoccupied band increasing faster 
than the gap in the initially occupied band. As a re- 
sult, although for the undoped system these gaps have 
very different orders of magnitude, beyond /i = fi c they 
achieve comparable values, being manifested in the par- 
tial gaps as functions of \i and T (Fig. 02(a)). Our results 
are independent of the character of the inter-band in- 
teraction: for an attractive (repulsive) interaction, the 
gaps on two bands have the same (opposite) signs, but 
the thermodynamic properties discussed here remain the 
same. We also discuss possible connections with STO 
interfaces, proposing tunneling experiments on both the 
LAO/STO heterostructures and STO thin films to shed 
light on the superconducting state of these heterostruc- 
tures. Finally, we discuss the importance of the dimen- 
sionality, arguing that most of the sharp features ob- 
served at \i = /i c in the two-dimensional model may be 
smoothed out by three-dimensional band dispersions. 

Our paper is organized as follows: in Section [III we for- 
mulate and solve our two-dimensional two-band model, 
discussing how T c and the gap functions change as the 
chemical potential is varied. In Section IIII( we calculate 
the LDOS corresponding to the two-band model, com- 
paring our results to the experiments of Binnig et al Q/. 
The effects of three-dimensionality are discussed in Sec- 
tion IIV| and our concluding remarks are presented in 
Section El 



II. TWO-DIMENSIONAL TWO-BAND MODEL 

A. Formulation of the model 

First-principle calculations @, [Io[ reveal that the band 
structure of bulk STO near the Fermi level is composed 



of the Ti t2 g orbitals d xz , d yz , and d xy . Due to the spin- 
orbit coupling, there is a splitting of 2ji c ~ 20 — 30 meV 
between the d xz j yz and d xy orbitals, with the latter hav- 
ing a larger onsite energy than the former. There is also 
another splitting between the d xz and d yz orbitals, which 
is however one order of magnitude smaller than \i c . Here, 
we will neglect this extra splitting and consider an ef- 
fective "heavy" degenerate d xz / yz electron-like band 1, 
whose bottom is located fi c + fi below the Fermi level, 
and a "light" electron-like d xy band 2, whose bottom is lo- 
cated 2/i c above the bottom of the band 1 (see Fig. [T^a)). 
Notice that in the LAO/STO interfaces (shown in Fig. 
Gib)), there is an orbital reconstruction that shifts the d xy 
band below the d xz / yz [18(. While this will change the 
possible orbital occupations of the electrons, the effective 
two-band model considered here is still representative of 
the system. 

Introducing the electronic creation operators c^ kcr , 
with band index a, momentum k, and spin a, we 
have the following interacting Hamiltonian in the pair- 
ing (particle-particle) channel: 

H = Yj ( £a ' k ~ ^ C l,ka C «,ka + (1) 
a,~k,a 

Vab,kk'cl M c\ _ k ^C 6 _ k /^C 6>k / t + h.c. 

a,6,k,k' 

where e a ,k gives the band dispersion s\ y k = k 2 /2mi — 
l^c,£2,k — k 2 /2rri2 + /i c > I 1 is the chemical potential, and 
V a b,\zk.' are the pairing interactions. Hereafter, we assume 
them to be momentum-independent, i.e. we consider s- 
wave SC states and for convenience choose the zero of 
energy to be at the midpoint between the two energy 
band minima. Furthermore, we consider that the pairing 
involves only states near the Fermi level, yielding the 
mean- field 2x2 system of gap equations: 




Here, W is the upper cutoff of the interaction, A$ are 
the SC gaps, pi are the density of states at the Fermi level 
(in 2D the DOS is assumed to be independent of energy), 
111 = —ii— ii c is the position of the bottom of the electron- 
band 1, and fi2 = li c — M is the position of the bottom 
of the second band. For simplicity, we will discuss our 
results in terms of the dimensionless coupling constants 
\j = —PiVij- We assume V21 = V12 and therefore A21 = 

(P2/Pl)Al2- 

From Eqs. ([2j) it is straightforward to obtain T c . One 
linearizes the gap equations and calculates their largest 
eigenvalue, which equals 1 at T c . This procedure yields 
the implicit equation: 
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Figure 2: (Color Online) (a-b) Superconducting transition 
temperature T c and (c-d) the T = ratio of the gap as a 
function of chemical potential for various values of A 12 and 
A22. The light gray dotted line refers to A2/A1 = 1. Note, 
the gap of the second band is finite at all doping, but it is 
significantly smaller for \i < fi c . In all panels, An = 0.14. 
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We note that T c <C /i c <C W, which justifies the use 
of the log functions above. The last term in I2 (//) is 
the one responsible for the changes in T c as the second 
band crosses the Fermi level. In particular, very close to 
the point \i = /i c , this function has a linear dependence 
on /i — /i c , whereas far from this point, it behaves as the 
logarithm of \i — fi c . We can also obtain the gap functions 
Ai from Eqs. (j2j) in a straightforward way. For T = 0, we 
obtain simple algebraic equations for A^o- For non-zero 
T < T c , one has to solve the self-consistent gap equations 
numerically. 

Multi-band gap equations such as ([2]) have been inves- 
tigated in the past in the context of different compounds, 
such as iron pnictides and MgB2. Here, our interest is an 



aspect of these gap equations that has not been as widely 
investigated |30l [31], namely, the evolution of the solu- 
tions of the gap equations as a function of the chemical 
potential - particularly near the point where the second 
band starts to be occupied, [i — \i c . 



B. Results for T c and Ai 

We first need to discuss the values of the parameters 
that describe the electron-doped STO compounds. First, 
we note that the two bands have different orbital content, 
i.e. band 1 is composed of the d xz / yz orbitals, whereas 
band 2 is composed of d xy orbitals. Since intra-orbital 
interactions are typically larger than inter-orbital inter- 
actions, we expect the intra-band pairing to be larger 
than the inter-band pairing, i.e. |A^| ^> |A^-|. As we 
will show below, this gives results in agreement with the 
experimental observations of Binnig et al. 

Under these conditions, SC only appears for attractive 
intra-band interactions, A^ > 0, regardless of the sign of 
A12, since only \\ 2 appears in the equation ([3]) for T c . 
The sign of the inter-band coupling A12, however, affects 
the structure of the SC state, as the eigenvector that 
diagonalizes the linearized form of Eqs. ([2]) is such that 
sign(AiA2) = sign(Ai2). Thus, for A12 > (attractive 
inter-band interaction V12 < 0), Ai and A2 have the 
same sign (the so-called s ++ state), whereas for A12 < 
(repulsive inter-band interaction V12 > 0), Ai and A2 
have opposite signs (the s +_ state) [30]. Since the SC free 
energy depends only on A12A1A2, the thermodynamic 
properties discussed in this paper do not depend on the 
relative sign between Ai and A2. Therefore, we set A12 > 
0, but it is possible that STO is an unconventional s-wave 
superconductor with structure. 

In Fig. [2j we present the enhancement of T c as well as 
the ratio between the zero-temperature gaps A20/A10 as 
a function of fi/fi c for several different values of A12 and 
A22, with An = 0.14 fixed. Since band 2 is lighter than 
band 1, we have p2 = (^2/^1)^1 = 0.4pi, in accordance 
with band structure calculations [8| . For the upper-cutoff 
W. we set W = 10/i c , and we consider An w A22 <C 1 
as well as A12 <C A^, as discussed above. As it is shown 
in the figure, even when the inter-band pairing is small, 
the occupation of the second band at fi = /i c is marked 
by a sharp increase in T c and by an enhancement in the 
zero-temperature gaps A10 and A20. While for \i < \i c 
the gap A10 is significantly larger than A20, for \i > 
He they become of similar magnitudes. The complete 
temperature behavior of the gap functions is shown in 
Fig. [31(a), which corresponds roughly to two BCS-like 
curves. 

If we considered a larger inter-band interaction A12, 
then we would have obtained a larger enhancement of T c 
at fi = /i c , but the gap function A20 would be much closer 
to A10 even for \i < \i c (shown in Fig. [2](a) and (c)). 
Meanwhile, the tunneling experiments of Binnig et al. 
[l| suggest instead that A20 is considerably smaller than 
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Figure 3: (Color Online) (a) Superconducting gap as a 
function of chemical potential and temperature, displaying 
the standard BCS temperature dependence, (b) Calculated 
LDOS for the two-band model as function of increasing chemi- 
cal potential and energy, (c)-(e) The chemical potential levels 
with respect to the two bands (left) and the calculated LDOS 
for the two-band model for p/p c = (c) 0.5, (d) 1.01, and (e) 
1.5 and for the set of parameters An =0.14, A22 = 0.13, and 
A12 = 0.02 



A 10 near p = p c . Conversely, if we considered a smaller 
intra-band interaction A22 < An, then the enhancement 
of T c at p = p c would be rather weak, in disagreement 
with the same data. Therefore, unless some fine-tuning 
is claimed, our results put important constraints on the 
character of the multi-band SC of electron-doped STO, 
as it appears to be dominated by attractive intra-band 
pairing interactions and weaker inter-band pairing. 

Note that this situation is the opposite of other multi- 
band superconductors, such as the iron pnictides. These 
materials have a much larger repulsive inter-band interac- 
tion A 2 2 ^> A11A22 enhanced by spin fluctuations, which 
gives rise to a sign-reversal s +_ state (H, (30[. Contrary 
to the iron pnictides, in STO the SC state is not near 
other electronic ordered states whose fluctuations can en- 
hance A12 . 



N(u) = --Im V[z/iG n (cc;,k) + z/ 2 G22(^,k)] 

7T < ^ 



(5) 



where denotes the degeneracy of the band. Within 
the mean field approximation each Green's function takes 
the BCS form Gjj = (o;+^(k))/[o; 2 -^ 2 (k)-A J 2 ] with the 
respective gaps from Eq. O By replacing the summation 
over momenta by an integration over energies, the LDOS 
can be calculated analytically, yielding 



N(uj) = Re 



E 



(cj + if] + Cj) pjjcj) + (u - Cj) Pj(-Cj) 
2c, 



(6) 

where pj are the density of states of each band, Cj 

= ^J(uj + irj) — A 2 , and r] provides a small width 

to simulate experimental measurements. For the two- 
dimensional model considered so far, pj are constants. 

Using the T = gaps obtained from the solution of 
Eqs. ([2]), we calculate N (uj). We observe two different 
behaviors, as shown in Fig. [3f b). For p < p c: N (uj) 
has a sharp peak near uo = A10, and the smaller gap A20 
shows no noticeable signature. This is due to the fact 
that this gap is not open at the Fermi level, but at an 
energy near the bottom of the second band p c — p^> A20. 
We note that the LDOS has a clear increase at = p c — p 
due to the additional contribution of the states coming 
from the second band. Only when p c — p ~ A20, we see 
a weak signature of the second gap at energies slightly 
larger than A10, where a clear peak is still observed. 

For p > p Cl after the second band crosses the Fermi 
level, the second gap gives rise to a smaller peak at en- 
ergies lower than the first peak, since A20 < Aio- As 
p increases, this peak becomes larger and moves closer 
to the peak coming from the first-band gap. Eventually 
these two peaks become difficult to distinguish, since the 
gaps become comparable in magnitude. As a result, the 
LDOS has a sharp peak preceded by a shoulder-like fea- 
ture (shown in Fig. Efc)-(e)). These general features 
of the LDOS spectrum are in good agreement with the 
tunneling data of Ref . [l| . 



IV. THREE-DIMENSIONAL CASE 



III. TUNNELING SPECTRUM 

The existence of two SC gaps has important implica- 
tions for several thermodynamic quantities. For instance, 
two-gap superconductivity can be inferred from the tem- 
perature dependence of the specific heat and penetration 
depth. Yet, the most direct probe for multi-gap super- 
conductivity is the local density of states (LDOS), which 
can be measured via tunneling spectroscopy. 

Within our two-band model, the total LDOS N (uo) is 
the sum of the contributions from each band: 



So far we have focused on the 2D case, which is most 
relevant to STO thin films, interfaces, or probes sensitive 
to surface states, such as the tunneling experiments by 
Binnig et al. [1]. In the bulk STO material, the 3D dis- 
persion gives rise to a density of states p (£) ex y^? which 
is expected to smoothen the sharp features observed in 
T c and A^o at p = p c . Nevertheless, as we will show in 
this section, the enhancement of T c and A^o at p = p c is 
still present. 

We continue to assume parabolic bands. Then, 
the density of states of each band is given by pi = 
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Figure 4: (Color Online) Calculations for the 3D two-band 
model as function of chemical potential for An = 0.14 and 
A 22 = 0.13. We present T c (a) and A 2 /Ai (b) for A12 = 
0.02, 0.04, and 0.08. The light gray dotted line in (b) across 
A2/A1 = 1 denotes the point at which the second gap crosses 
the first gap. The individual energy gaps (A*) for A 12 = 0.02 
are shown in (c). 



Pi,oV e + M ± Mc and the gap equations are changed ac- 
cordingly to: 



3=1 



2 A /e 2 + A? 



(7) 



ie + a| 



2T 



The solution of this system of equations is straightfor- 
ward. In Fig. 2J we show the detailed behavior of the 3D 
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Figure 5: (Color Online) Calculated LDOS for the weakly 
coupled 3D two-band model as a function of cj//jl c for \ij \i c 
= 0.5 (solid-black), 0.8 (dashed-red), 1.2 (dotted-blue), 1.5 
(dash-dotted-green). Here, we use the parameters An = 0.14, 
A22 = 0.13, A12 = 0.02 and k = 0.4. 



two-band model as function of the chemical potential for 
An = 0.14 and A22 = 0.13 (the same parameters as in the 
2D case). The evolution of T c does produce an "upturn" 
around ji = /i c , as shown in Fig. HKa). The ratio of the 
gaps, presented in Fig. HJb), displays a kink as the second 
band is populated. However, as shown in Fig. H^c), the 
individual energy gaps display a distinct increase. This 
is consistent with the data observed by Binnig et al [l|, 
where T c and the gaps clearly increase as function of the 
chemical potential. The T c data in Binnig et al [l[ shows 
a slight kink as the second band becomes populated. The 
3D model does not explain this slight kink and may in- 
dicate the presence of a 2D band structure from surface 
states. Therefore, we conclude that a two dimensional 
model is more accurate to describe the STO systems and 
may be the leading cause for the creation of a 2D electron 
gas at the interface of LAO and STO. 

For completeness, we also calculate the local density 
of states (LDOS) resulting from the three-dimensional, 
two-band model. The LDOS is also given by Eq. (j6j), but 
with the appropriate energy-dependent density of states 
pi. Figure [5] shows the results, with the LDOS plotted 
as a function of Lj/fJL c f° r p/Pc = 0.5, 0.8, 1.2, 1.5. The 
behavior is similar to the one observed by Binnig et al. 
in Ref. [l| , where the population of the second band pro- 
duces a small lower- frequency shoulder and an increase 
in the LDOS. 



V. CONCLUSIONS 

The observation of multiple gap structures in the tun- 
neling spectrum of electron-doped bulk STO [l| is a well 
known fact, yet it has been largely overlooked in the dis- 
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cussion on the origin of superconductivity in both doped 
STO and LAO/STO intefaces. Motivated by this ob- 
servation we proposed a multi-band model for supercon- 
ductivity in electron-doped STO. The d-band manifold 
is split into two d xy and d xz j yz bands with a splitting of 
around 10 meV before doping. We show that with in- 
creasing electronic occupation, the higher d-band(s) be- 
come occupied, engaging an additional source of electrons 
in the pairing. The onset of the occupation of the higher 
bands results in an increased DOS and pairing conden- 
sate gain, which in turn results in a sharp increase in 
T c and in the SC gap in the second band. Comparing 
to experiments, we find the primary source of pairing to 
be the intra-band coupling An ,22 with smaller inter-band 
Josephson coupling A12 <C An ,22- Depending on the sign 
of the inter-band coupling, the SC gaps on the two bands 
can have the same or opposite signs, resulting in either 
«s ++ or SC states. Although our results were mo- 
tived by bulk STO, the fact that T c in the LAO/STO 
interfaces is very close to the T c of electron-doped STO 
@ , and the observation that the maximum T c of the het- 
erostructure takes place when a second d-band becomes 
occupied [25|, HH, we suggest that the same multi-band 
effects may take place in the oxide interfaces. Of course, 
in the heterostructures, the breaking of inversion sym- 



metry at the surface may lead to additional contribu- 
tions to T c , not considered in our model (32|. To test our 
proposal, we suggest to use simple direct spectroscopy 
in the SC state of doped STO thin films and LAO/STO 
interfaces, such as scanning tunneling microscopy or pla- 
nar tunneling. Both probes would reveal the existence 
of a two gap structure if indeed superconductivity is a 
two-band phenomena in these systems. The observation 
of two band features would likely be possible only in a 
clean limit when scattering will not wash away features 
of the weaker gap. 
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